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Abstract 

Photo-acoustic Tomography (PAT) and Thermo-acoustic Tomography (TAT) are medical 
' (— i '- imaging modahties that combine the high contrast of radiative properties of tissues with the 

high resolution of ultrasound. In both modalities, a first step concerns the reconstruction of the 
radiation-induced source of ultrasound. Transient Elastography (TE) and Magnetic Resonance 
Elastography (MRE) combine the high elastic contrast of tissues with the high resolution of 
ultrasound and magnetic resonance, respectively. In both modalities, a first step concerns the 
reconstruction of the elastic displacement. 

The result of this first step, which is not considered in this paper, is the availability of 
internal functionals of the unknown tissue properties. All imaging modalities are recast as the 
\ reconstruction of parameters in elliptic equations from knowledge of solutions to such equations. 

This paper provides a characterization of the parameters that may or may not be reconstructed 
from such internal functionals. We provide explicit reconstruction procedures and indicate how 
stable they are with respect to errors in the available measurements. 

The modalities PAT, TAT, TE, and MRE allow us to reconstruct high-contrast optical 
and elastic properties of tissues with the high resolution of ultrasound or magnetic resonance 
imaging. They provide a means to reconstruct second-order tensors modeling tissue anisotropy 
as well as complex- valued coefficients modeling absorbing and dissipative effects. 

1 Introduction 

The optical and electrical properties of tissues often display a high contrast between healthy and 
non- healthy tissues [1, 2, 3, 4]. However, the reconstruction of such coefficients from boundary 
measurements is mathematically an ill-posed problem. As a result, modalities such as Optical To- 
mography and Electrical Impedance Tomography are typically low-resolution. In recent years, several 
methods have been developed to combine the high contrast of optical- and electrical- based modalities 
with the high resolution of other, often low-contrast, modalities such as ultrasound tomography or 
magnetic resonance imaging. Photo-acoustic tomography (PAT) and thermo-acoustic tomography 
(TAT) are such modalities combining high optical contrast with high ultrasound resolution. 
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The physical mechanism that allows the coupling between radiation and ultrasound in PAT and 
TAT is the photo-acoustic effect, which may be described as follows. In both PAT and TAT, pulsed 
radiation propagating through a domain of interest is partially absorbed. The thermal expansion 
resulting from the absorption of radiation is responsible for the emission of ultrasound. An array of 
transducers then records the emitted ultrasound at the boundary of the domain. This is the available 
information from which we aim to reconstruct the optical coefficients of the tissues. 

In PAT, radiation is typically near-infra-red light, while in TAT, radiation is typically low fre- 
quency (on the order of hundreds of MHz) electromagnetic radiation [5, 6, 7]. In both PAT and 
TAT, a first step of the reconstruction procedure consists of reconstructing the amount of emitted 
ultrasound from the available measurements. Mathematically, this is an inverse wave problem that 
aims to reconstruct the initial condition from boundary time-dependent measurements. We assume 
this first step done; see e.g.: [8, 9, 7, 10, 11] for the bio-engineering literature and [12, 6, 13, 14] for 
the mathematical literature. 

The maps reconstructed during the first step provide important qualitative information. However, 
they do not quantify the optical properties of the tissues. This paper is concerned with the second, 
quantitative, step of PAT and TAT, called QPAT and QTAT, respectively; see, e.g., [8, 15, 16, 17, 18] 
in the bio-medical literature and [19, 20, 21, 22, 23, 24] in the mathematical literature. 

From the analyses in, e.g., [19, 22], we know that not all of the unknown tissue parameters 
may be reconstructed in QPAT experiments. This paper follows recent results in [22] to obtain a 
precise, general, description of all that can be reconstructed from QPAT and QTAT data. We also 
describe the reconstructions in the elasticity-based imaging modalities Transient Elastography (TE) 
and Magnetic Resonance Elastography (MRE) when they are modeled by a scalar equation; see 
[25, 26, 27, 28, 29, 30, 31, 32]. 

The elliptic models describing radiation propagation and the internal functionals resulting from 
solving the first step in QPAT and QTAT are presented in section 2. The analysis of what may or 
may not be reconstructed from such information is carried out in section 3. The main mathematical 
properties that allows us to carry out the analysis is that the ratio of internal functionals is equal to 
the ratio of the elliptic solutions. This property is satisfied by QPAT, QTAT, TE, and MRE. What 
may be reconstructed from ratios is presented in section 3.1. The final step of the reconstruction 
depends on the modality of interest. QPAT is addressed in section 3.2, QTAT in section 3.3, and TE 
and MRE in section 3.4. In QPAT and QTAT, not all the coefficients can be uniquely reconstructed. 
We consider several scenarios that lead to unique reconstructions. We stress that anisotropic as well 
as complex-valued coefficients can be reconstructed for all modalities. The results are discussed in 
the concluding section 4. 

QPAT, QTAT, TE and MRE belong to a class of imaging modalities now often referred to as 
hybrid inverse problems or coupled-physics inverse problems. We refer the reader to e.g. [33, 34, 35, 6] 
for recent references on the rapidly evolving field of hybrid inverse problems. 

2 Quantitative PAT and TAT 

The first step of QPAT and QTAT concerns the reconstruction of the map of ultrasound generated 
by the propagating radiation. We assume this first step done; see [8, 12, 6, 13, 7, 10, 11] for an up-to- 
date list of references on this problem. The second step of QPAT and QTAT provides quantitative 
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reconstructions of the optical parameters from knowledge of the map obtained in the first step. The 
mathematical problem associated to this second step is presented in section 2.1 for QPAT and in 
section 2.2 for QTAT. 

2.1 Quantitative Photo- Acoustic Tomography 

In PAT, the propagation of radiation is modeled by the following diffusion equation 

—V ■ 7Vnj + (TUj = in X, Uj = fj on dX. (1) 

Here Uj is the spatial photon density and (7, cr) are the diffusion tensor and absorption coefficient, 
respectively. The boundary conditions fj for 1 < i < J indicate the different ways the domain is 
probed. Here, we use Dirichlet conditions to simplify. The reconstruction procedures do not depend 
on the choice of boundary conditions, which can thus be replaced by Neumann or more general Robin 
(mixed) conditions. X is a regular, open, bounded, domain in where n > 1 is spatial dimension. 
As radiation propagates through the domain, the resulting acoustic signal is characterized by 

Hj{x) = T{x)(T{x)uj{x). (2) 

Here, F is the Griineisen coefficient, quantifying the coupling between the absorbed radiation energy 
a{x)uj{x) and the amount of resulting acoustic signal. 

The general problem of QPAT is therefore to know what can be reconstructed in 

(7,^,r) (3) 

from knowledge of Hj{x) for I < j < J. See, e.g., [23, 8, 7] for additional details. 

2.2 Quantitative Thermo- Acoustic Tomography 

In TAT, the propagation of radiation is modeled by Maxwell's equations 

V X ^ = -dtB, V xH = J + dtD, (4) 

where {E, B, H, J, D) are the electric field, the magnetic field, the electrical displacement field, the 
electrical current density, and the magnetic flux density, respectively. For linear media, which we 
assume here, we have the relations J = aE, D = eE, and B = fiH, where (j{x) is the conductivity, 
e{x) the relative permittivity, and fi{x) the relative permeability. After eliminations, we find the 
equation for the electric field 

-V X /i" V xE = adtE + ed'^E. (5) 

Replacing the above system of equations by a scalar model, we obtain formally the scalar equation 
V ■ Vu = adtu + ed^u. Here m, may be seen as one component of the electric field in a simplified 
setting. The above derivation may be justified when V ■ i? = and /i is constant. See [21] for 
partial results for (5). With time-harmonic sources and solutions with frequency the equation for 
u becomes the Helmholtz equation 

V ■ jj'^Vu + (ew^ — iua)u = in X, u = f on dX. (6) 
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As for PAT, more general boundary conditions can also be considered. 

As radiation propagates through X, the emitted acoustic signal is described by 

H{x) = T{x)(j{x)\u\'^{x). (7) 

As above, the Griineisen coefficient F describes the coupling between the absorbed radiation energy 
a{x)\u\'^{x) and the resulting acoustic signal. 

Let us now assume that the medium is probed in K different ways fj for 1 < j < K. Solutions of 
the form are sufficient to generate information of the form Ujul- Indeed, using the polarization 
formula for the inner product on C, we obtain that 

UkU*j = ^(kj + Uk\^ - i\uj + iuk\^ - (1 - i){\uj\^ + \uk\'^)y (8) 

Therefore, using four physical experiments with boundary conditions fj, fk, fj + fk, fj + ifk, we can 
reconstruct T{x)a{x)\uj\'^, r{x)a{x)\uk\'^, as well as T{x)a{x)ukUj. We thus assume the availability 
of the internal functionals 

Hj{x) = T{x)(j{x)ujul, l<j<J. (9) 

We refer to the above information (Hj) as a J— dimensional data set even though the number of 
physical experiments K necessary to acquire Hj may be larger than J. The inverse problem of 
QTAT therefore consists in understanding what can be reconstructed in 

{fi-\x),e{x),a{x),T{x)), (10) 

from knowledge of {Hj)i<j<j. See e.g. [36, 21, 7] for additional information about this model. 

2.3 General setting for QPAT and QTAT 

QPAT and QTAT thus aim to reconstruct the coefficients (a, c, F) from knowledge of 

Hj{x) = T{x)c{x)uj{x) in the QPAT setting, (11) 

Hj{x) = r{x)'^c{x)ujul{x) in the QTAT setting, (12) 

where Uj is the solution to the elliptic equation 

V • aWuj + cuj = 0, in X, Uj = fj on dX, ^ < j < J- (13) 

Here, is the imaginary part of c. The coefficients (a, c) are possibly complex-valued. 

3 Reconstruction procedure 

The functionals in QPAT, QTAT, TE, and MRE share the property that 

Hk Uk 

We exploit this structure to obtain preliminary reconstructions in section 3.1. The final step of the 
reconstruction is modality-specific. QPAT is treated in section 3.2 while QTAT is handled in section 
3.3. The elasticity-based modalities TE and MRE are described in section 3.4. 



(14) 
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3.1 Reconstruction from solution ratios 

We assume that experiments are carried out for J G N boundary conditions /j, 1 < i < </, with 
resulting solutions Uj of (13). In the first step of the reconstruction, we assume the availability of 
ratios of solutions. We make our first assumption that 

ui does not vanish throughout the domain X, (15) 

and define 

vj = ^, l<j<J-l. (16) 

Since the functionals Hj are linear in the solutions Uj in QPAT and are bilinear in («*, Uj) in QTAT, 
we observe that Vj is known in the QPAT and QTAT settings. 
Some straightforward algebra shows that 

a:V'^^Vj + {V-a)Vvj = in X, a{x) := a{x)ul{x). (17) 

Here, we use the notation: a : V®^fj = Xlfc i=i '^fc'^i'fc^^-i^i — Tr(aV®^i'j). We present the recon- 
struction of a from knowledge of a sufficient number of solutions vj following [22]. 

Reconstructions in the presence of a scalar diffusion coefficient. Let us assume that a is 
scalar. Then (17) is equivalent to the equation 

Av^ + . Vv, = 0. (18) 

aui 

We make our second assumption: 

(Vf j)i<j<„ forms a basis of M"" at each point x E X. (19) 

We define Hij = Vvj ■ Vvj a matrix that is therefore invertible and the entries of the matrix 
H^^. For any complex- valued n-dimensional vector F, we then have the decomposition 

F = W^F -VvjVvi. (20) 

Here and below, we use the convention of summation over repeated indices [i and j above summed 
between 1 and n). Therefore, with (em)i<m<„ the standard orthonormal basis in M", we have 

em = (j)T'^Vi for <PT = W'em-Vv,\/v^. (21) 
This implies that for all 1 < m < n, we have 

V(a?/?) . ^ . Viauj) 

em ■ 2 + C^^J = 0' i-e-, 2 = -^rAt-j-em. 

auf auf 

In other words, the vector "^^^Jj^ is known explicitly. This is a redundant system of first-order 
equations for auf. For instance, choosing 1 < m < ra, we can solve the ordinary differential equation 

em -Viauj) + Av,){aul) = 0, 



provided that auf is known on dX. This provides an exphcit reconstruction procedure for auf. 

Reconstruction of a general, possibly anisotropic, diffusion tensor. The reconstruction 
when a is a full tensor requires a larger number of ratios Vj. We assume again that (Vwi, . . . V«„) is 
a basis of at each point x & X and define Hij = Vvi ■ Vfj with the entries of H~^. Let 

in(n + 3), (22) 
1<J<4-1, l<m<M^. (23) 

These expressions for the coefficients 9"^ ensure that T^^IT^ O^Vvj 
symmetric matrices 

/n-l 

M"^ = ^ ^rnV®2y., 1 < m < M„ 

We deduce from (17) that 

a:M'^ = Tr(aM'") = 0, 1 < m < M„, (25) 

for all X & X. We now make our third assumption: 

We assume that (M''")i<m<M„ forms a free family of symmetric matrices. (26) 

Since the dimension of the linear space of symmetric matrices is M„ + 1 = \n{n + 1), we deduce 
from (25) that a is in the orthogonal complement to the span of the M„ matrices M"*, which is a 
one-dimensional space. Let us call M°(x) a non-trivial matrix in that space. Such a matrix M° may 
be obtained by Gram-Schmidt orthogonalization [37] for instance observing that (/, {M'^)i<m<M„) 
forms a basis of the linear space of complex-valued, symmetric, matrices. This proves that 

a{x) = aulix) = r(x)M°(x), (27) 

for a scalar function r(x) to be determined. Using (20) and (17), we deduce that 

V ■ (tM°) = -rH'^f : V®\jVt;,. 

Since a, and hence M^, is invertible, we can recast the above equation as 

Vr + (M°)"^(V ■ M° + H'^M^ : V^^Vj\/Vi)T = 0. (28) 

This is a redundant system of linear first-order equations for r as in the case of a scalar. Knowledge 
of r at one point, for instance on dX, allows one to uniquely and stably reconstruct t on X. We 
have used In — I = n + Mn ratios {vj)i<j<i^^i to obtain a unique, explicit, and stable reconstruction 
of the tensor awf . 

Additional information does not provide new independent information. At this stage, we 
have reconstructed a = au\. Note that an^ = au\v'j_i is known as well for 2 < j < J. 



Mn = -n{n + !)-!, /„ = n + 1 + M„ 

r -W'Vvn,+n -Vvk l<3<n 
ef = I 1 j = n + m 

I otherwise 



= 0. Let us now construct the 

(24) 
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We now show that additional internal functional do not provide any new information when 
HjUk = HkUj. Indeed, let Uk correspond to a new boundary condition fx and let Hk be the 
corresponding internal functional. Then we find that 

V ■ aujV^ = in X, ^ = fy on dX. (29) 
^1 J^i Ji 

This is an elliptic equation on X with known Dirichlet conditions on dX. As a consequence, since 
a = au\ is known, then can uniquely be determined from the above equation. There is no need 
to acquire experimentally. 

A change of variables. All the QPAT, QTAT, TE or MRE information is thus encoded in {au{, Hi). 
We wish to recast au\ as a more explicit functional of the unknown coefficients (a, c). Let decompose 
a as a = B'^a with a a, possibly complex- valued, diffusion tensor such that det a = 1. The amplitude 
of a is written as 5^, where B could again be a complex- valued scalar. It is defined uniquely by 
continuity on the simply connected domain X. Since au\ is known and Ui is a scalar, then a is known 
as well. Let us now define v = Bui. Some algebra shows that 

V-aVB c 

Aw = qv on A, q = h (30) 

Note that = 5^ u\ is known since v'^a = auf, which is known. As a consequence, q is known. 

A summary. At this stage, we have used the equations for Uj and the information Vj for 1 < j < J—1 
to reconstruct aul. The number of necessary internal functionals is J = rz, + 1 in the case of a scalar 
coefficient a and J = In = \n{n + 3) in the case of a symmetric tensor a. 

All internal functionals Hk can be reconstructed from knowledge of Hi and au\., and therefore no 
additional information can be obtained by acquiring more measurements. We have then decomposed 
a = B'^a and used the equation for ui to eliminate it and derive knowledge of q in (30). All the 
QPAT, QTAT, TE, or MRE information available about the coefficients is encoded in {a,q,Hi). 
What we can extract from (a, q, Hi) now depends on the modality under consideration. 



3.2 Reconstructions in QPAT 

We recall that Hi{x) = T{x)c{x)ui{x). Knowledge of (a, q, Hi) is thus equivalent to that of 

(a,x,g)= a, — , + — . (31) 



5' B 

The above information (a, x, q) (i) has been reconstructed uniquely and stably from the available data; 
and (ii) is all that can be obtained about (a,c, F). This is the main result of this paper concerning 
QPAT. Note that two scalar functions, (x, q) are known while three scalar functions {B, c, T) are 
unknown. It is therefore impossible to reconstruct all of (a, c, F) from QPAT data without further 
prior assumptions. 

Following the derivation in [19], we obtain that knowledge of one function in {B,c,T) uniquely 
and stably determines the other two functions. For instance, if the Griineisen coefficient is assumed 
to be known, then (31) provides the following elliptic equation for B: 

V ■ aVB -qB + ^ = in X, (32) 
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with known boundary conditions on dX. As in[19], we can show that the above elhptic equation 
admits a unique solution. Alternatively, as also shown in [19], we find the equation 

V ■ aulV— + - = in X, — = — oxidX. (33) 

ui r ui fi 

Since auf is known, this is an elliptic equation for ui from which we then easily deduce a and c. 

We thus obtain that knowledge of F and QPAT data uniquely and explicitly determines the 
coefficients (a,c). Note that a is a possibly anisotropic tensor. When F is not known, then the 
anisotropy a can still be explicitly reconstructed. However {B, c, F) are reconstructed up to any 
transformation that leaves (XtQ) above invariant; see however [20]. 

3.3 Reconstructions in QTAT 

In QTAT, the coefficient c is naturally complex valued. In the general case, a is also possibly complex- 
valued while F is a positive, real- valued, coefficient. From Hi{x) = T {x)'^c{x)\ui\'^ (x) we obtain that 
the available information about (a, c, F) is 

X, q) = («, + — ) , (34) 

and this is all that can be reconstructed from QTAT data. Thus, the coefficients (a, c, F) can be 
reconstructed up to any transformation that leaves the above coefficients (a, x, q) invariant. 

When a is real-valued, then the imaginary part of q is '^cB~^. This means that the Griineisen 
coefficient F is then uniquely determined. The other coefficients {B, 3?c, Q^c) are determined only up 
to any transformation that leaves q invariant. 

In the QTAT setting, a = In practice, /i is always assumed to be a constant, known, scalar, 
parameter. In that setting, the available information is 

(FQ=c, 3?c, S3c) = (Fojcr, w^e, WO"). (35) 

We thus obtain that the three coefficients (T,e,a) can uniquely be reconstructed from J = n + 1 
QTAT functionals when is known. 

The stable reconstruction of a from J = 1 functional provided that (F, e) are known was proved 
in [21]; see also [36] for an explicit, less stable, formula from J = n + 1 functionals. 

3.4 Reconstructions in TE and MRE 

The mathematical tools presented above may be used to reconstruct anisotropic, complex-valued, 
coefficients in the imaging modalities Transient Elastography (TE) and Magnetic Resonance Tomog- 
raphy (MRE). As induced elastic waves propagate through the domain X, the resulting displacements 
are imaged either by ultra-fast ultrasound tomography in TE or by magnetic resonance imaging in 
MRE [25, 26, 27, 28, 29, 30, 31, 32]. Assuming a scalar model for the elastic displacement, m is a 
solution to the elliptic model (13) and we can consider that Hj{x) = Uj{x). 
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As a consequence, the reconstruction of aul yields that of a. Once a and ui are known, it is then 
straightforward to use (13) and get 

c(x) = (36) 

Ml 

In the TE and MRE settings, we thus observe that general complex-valued coefficients modeling 
possible dispersive effects (a, c) can be uniquely, explicitly, and stably reconstructed. Moreover, the 
coefficient a is allowed to be an arbitrary (elliptic) symmetric tensor. 



4 Discussion 

We have obtained in (31) and (34) a precise characterization of what can and cannot be reconstructed 
from a sufficiently large number of QPAT and QTAT experiments when radiation propagation is 
modeled by the scalar elliptic second-order equation (13). The main ingredients of the derivation are 
the structural property ^ = ^ of QPAT and QTAT functionals and the reconstruction procedure 
developed in [22]. Such a structure is also valid in TE and MRE, in which all the coefficients (within 
the approximation of a scalar model) are uniquely and explicitly reconstructed. 

Number of measurements. When the three assumptions (15)-(19)-(26) above are satisfied, J = 
n + 1 functionals are necessary to perform the reconstruction when the diffusion coefficient is scalar, 
whereas J = In = |n(n -|- 3) when a is a tensor. In fact, the reconstructions generalize [22] to the 
setting 

V ■ aVu + b-Vu + cu = x e X, (37) 

with b a complex-valued vector field. The above equation is not modified when (a, b, c) is replaced by 
(ra, rb — aVr, rc) for an arbitrary non-vanishing scalar function r. As a consequence, the number 
of degrees of freedom we can reconstruct in (a, b, c) equals |n(n-|-l)-|-n-|-l — 1 = /„. The number 
of internal functionals J = In is therefore natural in that setting and we obtain that J„ functionals 
precisely allow us to reconstruct /„ explicit functionals of the coefficients. 

The required number of internal functionals is therefore optimal when b ^ 0. However, the 
number of measurements in the general case and n -|- 1 in the case of a a scalar, may not be 
optimal. For instance, one coefficient in QTAT is reconstructed from one internal functional in [21] 
and two coefficients in QPAT are reconstructed from two internal functionals in [19, 23] independent 
of spatial dimension n. 

Constraints on the elliptic solutions and boundary conditions. The explicit reconstruction 
procedure requires that the three assumptions (15)-(19)-(26) be satisfied. These hypotheses have 
been shown to hold in [22] for "well-chosen" boundary conditions (/i, . . . , fj) in some specific situ- 
ations. In the generality considered in the present paper, the above properties are always satisfied 
locally [22]. The procedure presented in section 3 can then always be carried out locally on subsets 
of X. The reconstruction on the whole domain X may then require more functionals Hj than the 
number J; see [22] for additional details. 

Stability estimates. Hybrid inverse problems are being analyzed because they provide high reso- 
lution reconstructions. Stability estimates describe how errors in the acquisition of the functionals 
Hj propagate into errors in the reconstructed coefficient {a,x,Q)- 
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Stability estimates depend on the number of coefficients one wishes to reconstruct. For instance, 
in the QPAT setting with a and F real-valued and known and c = a real-valued, then a very simple 
procedure allows us to reconstruct a in a stable fashion. Indeed, let us assume that wi > cq > 
in the domain X. Then aui = ^ is known and hence Ui can be solved from the well posed 
elliptic problem (13). It remains to evaluate = ^ and we find the existence of a constant C 
such that ||5cr||oo < C'||5i/i||oo, where 5Hi is the error in the functional Hi and 5(t is the error in 
the reconstruction of the absorption coefficient. Here, errors are measured in the uniform norm 
ll/lloo = sup^jg^ 1/(3^)1 although the result also holds for other choices. 

In the QTAT setting, the reconstruction of a from one measurement i/i, provided that F and ^ 
and e are constants, has been analyzed in [21]. Again, we find that ||5cr||Y < C||5if||y for Y a space 
of sufficiently smooth functions; see [21]. The explicit iterative method in [21] is based on using a 
Banach fixed point. An explicit reconstruction from the n + 1 internal functionals Hj was recently 
presented in [36]. However, the reconstruction seems to involve a loss of three derivatives, whereas 
the methodology in [21] requires no such loss. With the same internal functionals, we reconstruct in 
this paper the three real- valued coefficients in (35) with a loss of one derivative (see below). 

In the QPAT setting with multiple unknown coefficients, the situation is less favorable. When a 
is not known, then (%, q) need to be reconstructed ffist. We refer to [19, 23] for different stability 
estimates for (x, q) when a is scalar and all coefficients are real- valued. When F is known, we find 
following [23, Theorem 4.1] that 

ll<^a||c"=(x) + < ||((5iJ)j||cfc+i(X;M2n) (38) 

for k > 2. Here ||/||c''=(x) = ^^Po<j<k,xex \f^''\^)\ with f^^^ the j—th. derivative of /. The recon- 
struction of {a, a) thus involves differentiating the data H = {Hi,H2) once. Such estimates still 
indicate that reconstructions should be accurate. The good behavior of the reconstruction of (a, cr) 
was conffimed by numerical simulations conducted in [19, 20]. 

Stability estimates for the general reconstructions presented in section 3 can be obtained following 
the derivation in [22]. The main conclusion of such estimates is that the reconstruction a full tensor a 
is less stable than when the coefficient a is scalar. Upon inspection of (24), we observe that the data 
Vj = -j^ need to be differentiated twice when is constructed. The reconstruction of Cm in (21) 

when a is scalar involves only first derivatives of the data Vj = (the estimated second derivatives 
are followed by one integration, which cancels the loss of one derivative). Consider the QPAT setting 
with F known. Then we find the stability estimates: 

\\Sa\\coix) + \\H\co(x) + ll^^llcMX) < C\\{SH)j\\c2(^x;RJ)- (39) 

We lose two derivatives to reconstruct (a, a). Since B is reconstructed using (28), we actually gain 
one derivative after integration and the stability estimates predict a better reconstruction of the 
scalar component B of a than the absorption coefficient a. 

Numerical implementation. The reconstructions presented above for QPAT when only (a, a) 
is unknown and a scalar and for QTAT when only a = 3?c is unknown have been implemented 
numerically in [19, 21]. In [19], the stability estimate (38) predicts that data need to be differentiated 
once. Of course, the differentiation of noisy data should not be done without some processing. A 
low-pass ffiter, for instance by a convolution with a kernel with adapted width, may for instance 
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be applied to the data prior to differentiation. The treatment of noisy data has been the object of 
considerable research in the inverse problems community. We refer the reader to e.g., [38, 39, 40] for 
a large class of standard methodologies to address the differentiation of noisy functions. 
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